I=imread('orion.jpg')
for i=1:900
for j=1:1200
smoothflt=ones(15,15)/225;
out=imfilter(I,smoothflt);
imshow(out)
for i=1:400
for j=1:296
if(out(i,j)>65)s(i,j)=255
else
s(i,j)=0;
end
end
end
imshow(out)
Out=I-out
imshow(Out)
OUT=Out+I;
imshow(OUT)